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This invention relates to conversion of images 
from gray scale pixel values to a reduced number of 
levels of pixel values, using error diffusion techniques 
that reduce visible artifacts normally noted with error 
diffusion. 5 

Image information, be it color or black and white, 
is commonly derived by scanning, initially at least, in 
a gray level format containing a large number of lev- 
els, e.g.: 256 levels for black and white and more than 
1 6 million (256 3 ) levels for color, which is usually un- 10 
printable by standard printers. In the following, the 
term "gray level" is used to described such data for 
both black and white and color applications. Standard 
printers print in a limited number of levels, either a 
spot or a no spot in the binary case, or a limited num- 15 
ber of levels associated with the spot, for example, 
four in the quaternary case. Gray level image data 
may be represented in a very large number of values. 
Accordingly, it is necessary to reduce gray level im- 
age data to a limited number of levels, so that it is 20 
printable. Besides gray level image information de- 
rived by scanning, certain processing techniques 
such as those described for example in European Pa- 
tent Application No. 91309640.0 produce gray level 
pixel values which require conversion to a limited set 25 
of legal, or correct output values. 

One standard method of converting gray level 
pixel image data to binary level pixel image data is 
through the use of dithering or halftoning processes. 
In such arrangements, over a given area, each gray 30 
level pixel within the area is compared to one of a set 
of preselected thresholds, comprising a matrix of 
threshold values or a halftone cell. The effect of such 
an arrangement is that, for an area where the image 
is gray, some of the thresholds within the matrix will 35 
be exceeded, while others are not. In the binary case, 
the portions of the matrix, or cell elements, in which 
the thresholds are exceeded are printed as black, 
while the remaining elements are allowed to remain 
white. The effect of the distribution of black and white 40 
over the cell is integrated by the human eye as gray. 
Dithering presents problems, however, in that the 
amount of gray within an original image is not main- 
tained exactly over an area, because the finite num- 
ber of elements inside each halftone cell only allows 45 
the reproduction of a finite number of gray levels. The 
error arising from the difference between the thresh- 
old value and the actual gray level value at any par- 
ticular cell is simply thrown away. This results in loss 
of image information. Dithering creates significant so 
image artifacts because it ignores this error com- 
pletely. A well known example is the banding or false 
contour artifact that can be seen in smooth image 
areas. Here, the image input gray level varies 
smoothly over an area while the halftoned image has 55 
to make a transition from one halftone dot (gray level) 
to another. This transition can clearly be seen as a 
band running through smooth image parts. 




170 A2 2 

Algorithms that convert gray images to binary or 
other number of level images attempting to preserve 
gray density exist, and include among them error dif- 
fusion, as taught, for example, in "An Adaptive Algo- 
rithm for Spatial Grayscale" by Floyd and Steinberg, 
Proceedings of the SID 17/2, 75-77 (1976) (herein- 
after, "Floyd and Steinberg"). Another, more elabo- 
rate method would be the error diffusion techniques 
of US-A-5,045,952. Additional modifications to the 
error diffusion algorithm taught by Floyd and Stein- 
berg have been proposed, e.g.: a different weighting 
matrix, as taught, for example, in "A Survey of Tech- 
niques for the Display of Continuous Tone Pictures on 
Bilevel Displays" by Jarvis et al., Computer Graphics 
and Image Processing, Vol. 5., pp. 13-40 (1976), and 
in "MECCA - A Multiple-Error Correction Computa- 
tion Algorithm for Bi-Level Image Hardcopy Repro- 
duction" by Stucki, IBM Res. Rep. R21060 (1981). 
Modifications of the error calculation and weight al- 
location have been taught, for example, in US-A- 
4,924,322 to Kurosawa et al., US-A-4,339,774 to 
Temple, and US-A-4,955,065,to Ulichney. 

Error diffusion attempts to maintain gray by mak- 
ing the conversion from gray pixels to binary or other 
level pixels on a pixel- by-pixel basis. The procedure 
examines each pixel with respect to a threshold, and 
the difference between the gray level pixel value and 
the threshold is then forwarded to a selected group 
of neighboring pixels, in accordance with a weighting 
scheme. The thus corrected image pixels are then 
considered input to the processing. In this way, the 
error calculated includes all errors previously made. 

An alternative error diffusion procedure is taught 
in "Images from computers" by M. R. Schroder (some- 
times spelled Schroeder) in IEEE Spectrum, pp 66-78 
(1 969) (hereinafter Schroder ). In this method the er- 
ror is only calculated between the original input pixel 
and the output, neglecting all previously made errors. 
This method leads to a poorer gray level representa- 
tion than Floyd and Steinberg but to higher image 
contrast. Modifications to the algorithm by Schroder 
are taught, for example, in "Design of Optimal Filters 
for Error-Feedback Quantization of Monochrome Pic- 
tures" by Jung Guk Kim and Gil Kim, Information Sci- 
ences 39, pp 285-298 (1986). 

In representing small variations in image appear- 
ance with an error diffusion technique, rather large 
areas of the image are required for the variation to be 
represented. Thus, while error diffusion is effective in 
maintaining gray over the image, it requires a rather 
large area of the image for the error compensation to 
be effective. Over such areas, undesirable textures 
and long range patterns, inherent in the distribution of 
error, may occur. 

It is believed that as the number of output levels 
increases, the need for error correction is reduced. 
Thus, error correction is essential in quantizing to a 
binary level, while less necessary for several level 
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systems. 

It is an object of the invention to provide a method 
and apparatus for the application of under compen- 
sated error diffusion to a pixel quantizing step in the 
conversion of image data from a number of input gray 5 
levels that is relatively large with respect to a number 
of desired output levels. 

According to the present invention there is pro- 
vided a method of quantizing pixel values in a black 
and white or color image, in which the image or each 10 
separation of a color image is formed by a plurality of 
pixels, each pixel representing an optical density of 
the image at a location within the image, and having 
an original optical density value associated therewith 
selected from one of a set of *c w original optical den- is 
sity values, where M c" is a number of members larger 
than a desired output set of "d" optical density values, 
characterised by: comparing the input optical density 
value of each pixel with the possible optical density 
values of the output and selecting one of these output 20 
values to represent the pixel; determining an error- 
term that is a difference between the new optical den- 
sity value and the original optical density value; deter- 
mining a preselected fraction of the error term to be 
used as the quantity for error correction; applying 25 
said preselected fraction of the error term on a 
weighted basis to each optical density value in a pre- 
determined set of neighboring pixels, thereby creat- 
ing modified optical density values of said pixels; and 
using said modified optical density values as input 30 
optical density values. 

Thus, gray pixel values, where the value of the 
pixel has one of a relatively large number of input lev- 
els, are directed through an error diffusion quantiza- 
tion process, for the conversion to one of a relatively 35 
small number of values. Each pixel of the gray level 
data is modified by a correspondingly weighted error 
correction term or terms from previously processed 
pixels, generating a modified pixel value. This modi- 
fied pixel value is compared to a set of threshold val- 40 
ues, where the comparison determines that the out- 
put will be one of a limited number of output levels. 
The difference value between the modified pixel val- 
ue and the output pixel value is distributed in accor- 
dance with a weighted distribution scheme through a 45 
matrix comprising a set of neighboring unprocessed 
pixels, increasing or decreasing the gray values of 
such pixels, which are referred to as modified pixel 
values. However, only a predetermined portion of the 
difference error is distributed, so that the effect of the so . 
quantizing error is limited. Image appearance is im- 
proved by limiting the extent of texture variation and 
long range patterns commonly associated with error 
diffusion techniques. 

In accordance with another aspect of the inven- 55 
tion, error may be transferred in the described man- 
ner either one dimensionally, to the next pixel in the 
scan line, or two dimensionally, to a preselected 



group of pixels including at least one other scan line. 

In accordance with yet another aspect of the in- 
vention, transfer of about half the quantization error 
from a single pixel has proven to have aesthetically 
desirable results for a 5 output levels system. How- 
ever, the amount of error to be transferred for aes- 
thetic improvement is related to the number of output 
levels, so that an increasing amount of error distribu- 
tion is desirable for systems with very few levels of 
output, while less error distribution is desirable for 
systems with larger numbers of output levels. 

In accordance with still another aspect of the in- 
vention, there is provided a method of quantizing 
gray level pixels, each pixel \ n j defined at one of 256 
levels of optical density, where the term 'optical den- 
sity* will be used hereinafter to describe the set of val- 
ues comprising the image data, be it density, intensi- 
ty, lightness or any other image data representation. 
Each pixel of the gray level data is modified by a cor- 
respondingly weighted error correction term or terms 
e„ f from previously processed pixels, generating a 
modified pixel value. This modified pixel value l n ,f + 
e nt f is quantized to one of five possible output levels 
using a predetermined distance measure. Usually 
this is done by comparing it to a set of four threshold 
values, where the comparison determines that the 
output value B n / representing the pixel in considera- 
tion will be selected from one of five output levels. The 
difference value e m between the modified pixel value 
l nf r + £n if and the output pixel value B n , f is distributed 
in accordance with a weighted distribution scheme 
through a matrix comprising a set of neighboring un- 
processed pixels, increasing or decreasing the gray 
values of such pixels, which are referred to as modi- 
fied pixel values. However, only a predetermined por- 
tion of the error e m is distributed, so that the effect of 
the quantizing error is limited. In accordance with one 
further aspect of the invention, a first effective error 
distribution scheme distributes half the quantizing er- 
ror e m to adjacent pixel (n + 1,?). In a second effective 
error distribution scheme half the quantizing error is 
distributed to a set of pixels in the adjacent scanline, 
and for pixel l n ,f, error light be distributed in accor- 
dance with 0.125e m to pixel (n + 1, f + 1), 0.25s m to 
pixel (n, t + 1), and 0.125e m to pixel (n-1, f + 1). 

These and other aspects of the invention will be- 
come apparent from the following descriptions to il- 
lustrate a preferred embodiment of the invention read 
in conjunction with the accompanying drawings in 
which: 

Figure 1 represents a prior art matrix for error dif- 
fusion weighted distribution; 
Figure 2 is a block diagram of a basic system for 
carrying out one embodiment of the present in- 
vention; and 

Figures 3A and 3B represents examples of two 
error diffusion matrices appropriate for use to- 
gether in the reduction of gray levels of image, in 
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accordance with the invention. 

Referring now to the drawings where the show- 
ings are for the purpose of describing the preferred 
embodiment of the invention and not for limiting 
same, a basic system for carrying out the present in- 
vention is shown in Figure 2. In the present case gray 
level image data from image input 8 is characterized 
as image data, each pixel of which is defined at a sin- 
gle level in a set of levels, the number of members in 
the set of levels being larger than desired. Each pixel 
will be processed in the manner described hereinbe- 
low, to redefine each pixel in terms of a new, smaller 
set of levels. Here, color data might be represented by 
a number of independent channels which are handled 
independently, or the color data might be represented 
as vectorial data in a predefined color space, e.g.: 
RGB. CIELab etc.. being submitted to vector opera- 
tions in the thresholding, error calculation and cor- 
rection. One common case of this method includes 
the conversion of data from a relatively large set of 
values to a smaller set of values for printing or other 
display technique. 

An input image of the type to be processed as 
hereinafter described, may be represented by a ser- 
ies of gray values (gray level pixels) arranged in an ar- 
ray of L hues, each hue containing N gray value pixels. 
As used herein. I n refers to a particular pixel in an im- 
age at position n. and the gray value or intensity lev- 
el of that pixel. The gray values are typically ex- 
pressed as integers with one example falling in the 
range from 0 to 255, although greater or lesser num- 
ber of levels are possible. An output image is consid- 
ered to consist of pixels, each pixel corresponding to 
an output dot or element that is printed by a digital 
printer. In one embodiment that will be herein consid- 
ered, image data is derived at 256 levels, and is quan- 
tized using the inventive method to one of 5 levels. 
Then, pulse width modulation is used to convert the 
5 levels of image data to a binary output, suitable for 
printing with an electrostatographic or xerographic 
laser printing device having binary output In electro- 
statographic or xerographic printing devices, the sin- 
gle small dots usually produced using error diffusion 
are not well reproduced. Using the pulse width mod- 
ulation method to convert the 5 levels of image data 
to a binary output allows pixels to be clustered togeth- 
er for better print output. Other methods of quantizing 
the 5 levels of image data to binary levels, which take 
into account the values of adjacent pixels, including 
general dithering, may be used to produce a set of e 
printer output values, where, in one preferred embodi- 
ment the number of members in the set is 2. Of 
course, the present invention has value independent 
of such an additional quantization step, particularly in 
display technologies, and in printing technologies 
such as ink jet printing, where the print dot can be well 
reproduced. 

With reference to Figure 2, the image input 8, 



which may be any source of gray level image data, in- 
troduces input image I into the system on a pixel-by- 
pixel basis. Each input pixel has a corresponding er- 
ror value e n ,f added to the input value l nf at adder 1 0, 

5 where e n ,f is the sum of error values of previous pix- 
els, resulting in a modified image, represented by 
modified pixel values, temporarily stored at block 12. 
The modified image, the sum of the input value and 
the error value of previous pixels (l n .r + £n ( r). is passed 

10 to threshold comparator 14. The modified image data 
is compared to threshold value(s) T, to determine an 
appropriate output level B n ,f for pixel \ n j and in the 
present embodiment, which of five levels B n ,f will rep- 
resent the gray level value. The threshold value(s) 

15 applied may be constant through the image or may 
vary either randomly, or in accordance with a dither 
pattern, as taught in "On the Error Diffusion Techni- 
que for Electronic Halftoning" by Billotet- Hoffmann 
and Bryngdahl, Proceedings of the SID, Vol. 24/3, 

20 (1 983), pp. 253-258, or in accordance with input infor- 
mation as taught in U.S. Patent Application 
07/396,272 entitled "Edge Enhanced Error Diffusion 
Algorithm u by R. Eschbach or in accordance with any 
other predetermined variation. Once output level B nt f 

25 is determined and directed to output image storage 16 
for eventual transfer to an output along line 18, the 
value of B nf is subtracted from the modified image 
value (ln f f + en,f) to generate the quantization error lev- 
el e m from quantizing pixel (l n ,f + e n .f). The subtraction 

30 operation is represented by the change sign block 20 
and subsequent adder 22, with e m representing the 
difference between the modified image value 
('n t f + en.f) and the output value B n f for pixel l n f . Error 
term e m is then distributed at error determination 

35 block 24 to a predetermined set of neighboring pixels 
in accordance with the scheme, that the given coef- 
ficients Ki-.-Knl 

0 < KiS'm + ... + K n e' m < e m 
where . . .K n are coefficients corresponding to the 

40 set of pixels to which error will be distributed, for the 
case of e m larger than zero. For e m smaller than zero, 
the "<" symbol is changed to a " > " symbol. 

There are several possible implementations of 
this scheme, two of which will be described here. In 

45 the first, the error value £,„ is multiplied by a constant 
a < 1 to generate e' m , using 

e'm = ae m and S,Kj = 1, 
the second implementation maintains Sm but uses 
Z,K,<1. 

50 The distribution matrix buffer 26 represents the ma- 
trix of stored values which will be used as weighting 
coefficients to generate e nt( as any pixel n, f that pass- 
es. 

In two example cases, in accordance with the in- 
55 vention, the matrices of Figures 3A and 3B show two 
error distribution arrangements. In Figure 3A, quan- 
tizing error is distributed one-dimensionally to the ad- 
jacent unprocessed pixel in the scan line in accor- 
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dance with a weighting scheme in which only 0.5e m is 
directed to the adjacent pixel. Considering that pixel 
l nt f+8 n ,f deviates from one of the possible output levels 
by not more than 12.5%. the error that can be passed 
is (0< I s m | <12.5%) x 0.5, consequently, the error re- 5 
ceived by pixel l n + 1nf is not more than ± 6.25% of the 
total dynamic range of the image, since all the error 
that is received by this pixel originated at pixel (n, f). 
The effect of this is to change the level of the next pix- 
el only when that pixel has a value that is very close 10 
to a threshold level or decision boundary. In accor- 
dance with this effect, the distribution of quantizing 
error has a limited cascading effect, beyond the initial 
pixel to which it is directed. Note that the numbers giv- 
en are for clarity only. The actual deviations depend 15 
on the number of levels and the quantization deci- 
sion. 

In Figure 3B, quantizing error is passed to three 
neighbors, as in the arrangement of Figure 2. There- 
fore for pixel l n P , error might be distributed in accor- 20 
dance with 0.125s m to pixel (n + f,f+ 7),.25e m to pixel 
(n,f+ f), and 0.125s m to pixel (n-f,f+ 1). Considering 
that pixel \ n j + 8 n ,f deviates from one of the possible 
output levels by not more than 12.5 % the maximum 
error that can be passed to pixel l„ + i, f + 1. is 25 
(0<E m <12.5%) x 0.125. Any pixel, however, receives 
error input from three different pixels, i.e. s n ,r the error 
to l n r is the sum of errorfrom !„_ f . l n f . and l n 
It can be seen that the maximum error received by l n f 
is again limited by 0.5 x (± 12.5%.) The effect of this 30 
is to change the level of the next pixel only when that 
pixel has a value that is very close to a threshold level 
or decision boundary. In accordance with this effect, 
the distribution of quantizing error has a limited cas- 
cading effect, beyond the initial pixel set to which it is 35 
directed. 

It will no doubt be appreciated that there are an 
enormous number of possible error distribution ma- 
trices available in error diffusion processes. What is 
required by this invention is that a predetermined por- 40 
tion of the error, less than 100%, is passed to the er- 
ror receiving set of neighboring pixels, for the reason 
of limiting the area over which an error compensation 
is attempted. 

It will also be appreciated that while the descrip- 45 
tion has referred to five level quantization, the sys- 
tems which provide multi-level output pixels other 
than five may use the invention as described with ap- 
propriate modifications to the total error feedback to 
accountforthe numberof pixel levels that can be pro- so 
duced in the output. Additionally, the system can be 
used to convert color data from a large set of possible 
color values to a smaller set of color values. 

The invention has been described with reference 
to a particular embodiment. Modifications and alter- 55 
ations will occur to others upon reading and under- 
standing this specification taken together with the 
drawings- The embodiments are but examples, and 



various alternatives, modifications, variations or im- 
provements may be made by those skilled in the art 
from this teaching which are intended to be encom- 
passed by the following claims. 



Claims 

1 . A method of quantizing pixel values in a black and 
white or color image, in which the image or each 
separation of a color image is formed by a plur- 
ality of pixels, each pixel representing an optical 
density of the image at a location within the im- 
age, and having an original optical density value 
associated therewith selected from one of a set 
of "c" original optical density values, where "c" is 
a number of members larger than a desired out- 
put set of "d* optical density values, characterised 
by: 

comparing the input optical density value 
of each pixel with the possible optical density val- 
ues of the output and selecting one of these out- 
put values to represent the pixel; 

determining an error term that is a differ- 
ence between the new optical density value and 
the original optical density value; 

determining a preselected fraction of the 
error term to be used as the quantity for error cor- 
rection; 

applying said preselected fraction of the 
error term on a weighted basis to each optical 
density value in a predetermined set of neighbor- 
ing pixels, thereby creating modified optical den- 
sity values of said pixels; and 

using said modified optical density values 
as input optical density values. 

2. The method of claim 1, wherein the predeter- 
mined set of neighboring pixels of pixel n, f, 
where n is the position of the pixel in scan line t, 
consists of the single pixel n + 

3. The method as described in claim 1, wherein the 
predetermined set of neighboring pixels of pixel 
n, f , where n is the position of the pixel in scan 
line ?, includes the set of pixels n.f + 1; n + 1,f + 
1;and n-1,f + 1. 

4. The method of claim 2, where the error term is 
applied in accordance with the weighting values: 
0.25 of the error term to pixel (n,f + 1 ), 0.125 of 
the error term to pixel (n + 1,f + 1 ), and 0. 125 
of the error term to pixel (n- 1,f + 1). 

5. The method of claim 1, where the preselected 
fraction of the error term to be applied to the pre- 
determined set of neighboring pixels is one half. 



BNSOOCID: <EP 0531170A2J_> 



i 



EP0 531 170 A2 



6. An apparatus for quantizing pixel values in an im- 
age formed by a plurality of pixels, each pixel rep- 
resenting an optical density of the image at a dis- 
crete location within the image, and having an 
original optical density value associated there- 5 
with selected from one of a set of "c" original opt- 
ical density values that has a number of members 
larger than a desired output set of "d" optical den- 
sity values, characterised by: 

a source of image data including pixels 10 
each having an original optical density value that 
is one of w c w original optical density values; 

an adder for determining the sum of a frac- 
tional error term derived from the quantization of 
at least one previous pixel and the original optical 15 
density value; 

a comparator for comparing the sum to at 
least one reference value, and producing an out- 
put value that is a member of the desired output 
set of d optical density values; 20 

an error term comparator for comparing 
the output value and the sum, and producing an 
error term that is the difference therebetween; 

means for determining a preselected frac- 
tion of the error term to be applied to each of a 25 
predetermined set of neighboring pixels; 

a memory storing the fractional error term 
to be applied to each original optical density value 
in a predetermined set of neighboring pixels; and 

means for directing the value of said pre- 30 
selected fraction of the fractional error term on a 
weighted basis to the adder for addition to each 
optical density value in the predetermined set of 
neighboring pixels. 

35 

7. The apparatus of claim 6, wherein the predeter- 
mined set of neighboring pixels of pixel n, P , 
where n is the position of the pixel in scan line £, 
includes the single pixel n + 1,f. 

40 

8. The apparatus of claim 6, wherein the predeter- 
mined set of neighboring pixels of pixel n, f, where 
n is the position of the pixel in scan line 2 , includes 
the set of pixels n.£ + 1, n + 1,f + 1, and n-1,P +1. 

45 

9. The apparatus of claim 8, where the error term is 
directed on a proportional basis to the adder for 
addition to the each optical density value in the 
predetermined set of neighboring pixels in accor- 
dance with the values: 0.25 of the error term to 50 
pixel (n,f + 1 ), 0.125 of the error term to pixel (n 

+ 1 J + 1), and 0. 125 of the error term to pixel (n- 

10. The apparatus of claim 6, where the predeter- 55 
mined fraction of the error term to be applied to 

the predetermined set of neighboring pixels is 
one half. 
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